Electron states, phonon-assisted relaxation and tunneling in self-assembled quantum 

dot molecules in an electric field 
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We present a theoretical analysis of the phonon-assisted relaxation in a system composed of 
two self-assembled vertically stacked quantum dots. We construct realistic model, which takes into 
account the geometry and strain distribution in the system. We calculate phonon-assisted relaxation 
rates between the two lowest states (in one- and two-electron cases). The relaxation rates and energy 
levels are studied as a function of external (axial) electric field and geometry of the structure (dot 
sizes). We show that the relaxation times can be as low as 1 ps but efficent relaxation occurs only 
for very finely tuned dots. 
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PACS numbers: 73.21.La, 73.63. Kv, 63.20.kd 



I. INTRODUCTION 

Coupled quantum dots are not only interesting from 
the theoretical point of view but also promising for de- 
signing quantum-coherent devices, including spin-based 
quantum bits. However, phonon-related processes are in- 
evitable in a cristal environment and may limit the feasi- 
bility of implementing quantum control in these systems. 
On the other hand, in many applications (e.g., lasers) 
fast relaxation between two states is desired. Therefore, 
it is essential to understand not only the spectral prop- 
erties of various kinds of semiconductor nanostructures 
but also phonon-induced processes. 

Carrier transfer and relaxation in self-assembled struc- 
tures was observed in many experiments. In those works, 
optical spectroscopy methods were used in stacked QDMs 
[THIS]. Miscellaneous mechanisms could be responsible 
for the observed properties. The heterogeneity of con- 
sidered structures and posible transfer mechanisms is re- 
flected in a relatively wide distribution of the measured 
transfer rates. However, in most cases the kinetics is at- 
tributed to tunneling [H [2 \7M El USHIS] ■ 

In recent works (Refs. [Ml [17]), we studied phonon- 
induced relaxation in one- and two-particle systems. In 
Ref. '^Tl, we investigated the dependence of the relax- 
ation on the system geometry. In this contribution, we 
also take into account dependence on external (axial) 
electric field for one- and two-particle systems. We show 
that phonon-induced relaxationmay be controlled either 
by system geometry or by external electric field. 



II. MODEL 

The system under consideration is composed of 
two vertically stacked QDs, formed in the Stransky- 
Krastanov self-assembly process. The dots are mod- 
eled as spherical domes of hights Hi and H2 (the index 
"1" refers to the lower dot) and base radii ri and r2- 
The domes rest on wetting layer of thickness ffwL and 
are separated by a distance D (base-to-base). The de- 
tails of the model and the parameter values are given in 



Ref. [T7]. We consider InAs dots in a GaAs enviroment. 
Approximately, our system has an axial symetry, so it 
is treated as cylindrically symmetric. This assumption 
strongly simplifies our calculations. The lattice constants 
mismatch leads to the appearance of strain. The strain 
distribution is represented by displacement fields that are 
found by minimizing the elastic energy of the system |16) . 
States of an electron in the strained nanostructure in 
the present approach are determined by the position- 
dependent conduction band edge and by the effective 
masses, which also vary across the structure. The lo- 
cal band structure is obtained from the 8-band Hamilto- 
nian with strain-induced terms (Bir-Pikus Hamiltonian) 
using the Lowdin elimination |161ll7j . From this, we cal- 
culate the conduction band edge and the effective mass 
tensor as a function of the spatial position. The electron 
wave functions are calculated within a variational multi- 
component envelope function approximation for equation 
with the Hamiltonian [16j 
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where m±{p, z),mz{p, z) are the components of the effec- 
tive mass tensor, E^i is electrical field magnitude, Ec{p, z) 
is the conduction band edge and e is electron charge. 

Next, we derive two-particle states within the standard 
configuration-interaction approach in a restricted basis 
of lowest-energy configurations (only two lowest single- 
electron states are included). The Hamiltonian of the 
interacting two-electron system has the form |il7) 
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FIG. 1. (a) The single-electron energy levels with a fixed size 
of both of dots (base radius r\ = r-^, = 10 nm) as a function 
of the electric field, (b) The single-electron energy levels with 
a fixed size of the lower dot as a function of the size (base 
radius r2) of the upper one. Here r\ = 10 nm, Ii\ — 3.7 nm, 
H2/r2 ~ 0.37. The energy reference level is 0.8 eV above the 
conduction band edge of unstrained bulk InAs. In both cases, 
the dot separation is D = 12 nm. 



Here fpiir) are single-particle eigenfunctions, £o is the 
vacuum permittivity, and £r is the dielectric constant of 
the semiconductor. After diagonalizing the Hamiltonian, 
we obtain three lowest singlet states. 

Finally, we derive the phonon-assisted relaxation be- 
tween the single-electron states as well as between the 
two lowest two-electron states. The coupling between the 
electrons and phonons is described by the Hamiltonian 
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where Fs^nn'io) is the coupling constants. We use the 
Fermi golden rule to calculate the relaxation rates be- 
tween the lowest energy eigenstates. In the calculation 
of phonon-induced relaxation rates, we take into account 
carrier-phonon couplings via deformation potential and 
piezoelectric interactions. 



III. RESULTS 

The single particle energy levels obtained within our 
approach for the fixed distance D = 120 nm between the 
dots are shown in Fig. [l|a-b). In Fig. [Tk, the energy is 
shown as a function of the external electrical field. Here, 
both dots have the same size. On the other hand, in 
Fig. [T]d, we show energy dependence on the upper dot 
size. The size of the lower dot is kept constant, while 
the base radius r2 of the upper dot and its height H2 are 
varied, with H2/r2 constant. Electronic (tunnel) cou- 
pling between the dots leads to the appearance of an 
anticrossing pattern. In Fig. [lb it occurs near the point 
where the dots become equal. In Fig. [T^ the resonance 
is located near the point Ec\ = 0. 

The spectrum of a two-electron system is dominated 
by antricrossings of different occupation configurations 



FIG. 2. The two-electron energy levels as a function of the 
electrical field (a) and of the upper dot size (b). The system 
parameters used in the calculations are identical as in the one- 
electron case. The energy reference level is 1.6 eV above the 
conduction band edge of unstrained bulk InAs. 
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FIG. 3. The thermalization rate as a function of the electrical 
field (a) and of the upper dot size (b). Results for a one 
electron system in a structure with r\ = 10 nm. Hi = 3.7 nm 
and H2/r2 = 0.37 at T = K (red solid line), 20 K (blue 
dashed line), and 40 K (green dotted line). 



PEl [TO] . In Fig. [2] the three lowest spin-singlet eigen- 
states of the interacting two-electron system are shown. 
In a similar way as in the one electron case, we present 
the dependence on the external electrical field (Fig. ^p) 
and on the upper dot size (Fig. [2k). The central res- 
onance occurs when the doubly occupied configurations 
have similar energy. This anticrossing is very narrow (less 
than 0.1 meV). The other two anticrossings occur at the 
degeneracy point between the singly occupied configura- 
tion and configurations with two electrons in the larger 
dot. These two splittings are wider than those appear- 
ing between the single-electron states, shown in Fig. [T] 
(2 meV vs. 1.5 meV). It is because the anticrossing 
of two-electron configurations is enhanced by Coulomb 
terms [j20 . 

The single particle relaxation rates are shown in Fig. [3] 
as functions of the upper dot size (Fig. [s}?) and the elec- 
trical field (Fig. [3^). These results are calculated for 
the same sample geometries as in Fig. [I] at three differ- 
ent temperatures. The overall magnitude of the relax- 
ation rates depends on the spatial overlap between the 



wave functions corresponding to the states involved in 
the transition. This is why the rates are large at the res- 
onance and become smaller when the system is shifted 
from the resonance point. Figs. |4^,b show the relaxation 



difference is that in the single electron case it tunnels to- 
wards an empty QD, whereas in the two electron case, 
there is already another electron. It leads to shifts (due 
to Coulomb interaction) to the parameter regimes where 
energies compensates for the Coulomb repulsion. 
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FIG. 4. The thermalization rate between the two lowest 
two-electron states at T = K (red solid line), 20 K (blue 
dashed line), and 40 K (green dotted line). Plots show the 
thermaUzation as a function of the electrical field (a) and of 
the upper dot size (b). 



IV. CONCLUSION 

We have studied the phonon-assisted relaxation (ther- 
malization) for single-electron and two-electron configu- 
rations in self-assembled quantum dots. Our model has 
taken into account the geometry and strain distribution 
in a QDM. The results show that the phonon relaxation 
is very efficient in a very narrow range of parameters 
near the anticrossings of energy levels. We conclude that 
relaxation times on the order of at least hundreds of pi- 
coseconds should be typical. The results for an electric 
field and various upper dot sizes are quantitatively consis- 
tent. Therefore, it is possible to control phonon-induced 
relaxation by system geometry as well as by an external 
electric field. 



rates for a two electron system. We use the same sam- 
ple geometries as in the one-electron case and investigate 
the dependence on electrical field and upper dot size. In 
general, the relaxation rates are similar to those found 
in the single electron case. This is due to the fact that 
both these processes are physically very similar. The only 
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